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AIRFOIL CHARACTERISTICS IN A NONEQUILIBRIUM AMBIENT 


STREAM BASED ON LINEARIZED THEORY 

By James J. Der 

Ames Research Center 
Moffett Field, Calif . 

SUMMARY 


The aerodynamic characteristics of a lifting body in a nonequilibrium 
stream are studied for a simple case, namely, the case of a nonequilibrium 
oncoming stream past a thin airfoil at a small angle of attack. Analytical 
solutions are obtained by the use of nonequilibrium linearized theory. 

Results show that when the relaxation length on the windward side of the air- 
foil is not appreciably different from that on the leeward side, a nonequi- 
librium ambient condition will normally induce a higher pressure and lower 
drag, but will not affect the lift, pitching moment, and the center of pres- 
sure. The local nonequilibrium condition behind the leading-edge shock, on 
the other hand, will affect all of the aerodynamic characteristics in general. 
When the relaxation length on the windward side is significantly different 
from that on the leeward side, on the other hand, the nonequilibrium ambient 
condition will also affect the lift, pitching moment, and the center of pres- 
sure. In general, the higher suction on the leeward side gives higher lift 
and larger pitching moment because the condition on the leeward side is 
closer to frozen. 


INTRODUCTION 


The present study is motivated by the desire to learn something about 
the effects of a nonequilibrium condition in the ambient stream on the aero- 
dynamic characteristics of a body. For an arbitrary body, such a study would 
normally require detailed numerical calculation of the flow around the body, 
with the free- stream condition out of equilibrium. Calculations of this type 
are, in general, tedious. 

A case for which we can obtain a simple analytical solution is the non- 
equilibrium flow past a thin airfoil at a small angle of attack. From the 
results of this case we can discover some indication of what the ambient non- 
equilibrium would do to the characteristics of a lifting body, with or without 
the local flow condition in equilibrium. 

In order to study problems whose solutions can be obtained analytically, 
we shall use a small- disturbance theory for nonequilibrium flow. This theory 
and its limitations have been described by Vincenti in reference 1, where he 
introduced a model for flow with nonequilibrium ambient conditions. The use 



of a small- disturbance theory limits the quantitative accuracy of the results. 
The results are not, therefore, intended to be accurate solutions to practical 
problems, but rather as qualitative information about the aerodynamic charac- 
teristics of airfoils in a nonequilibrium ambient stream. 

In the following sections, the physical aspects of the problem are first 
described. A solution to the governing equation is then obtained. The result 
is applied to biconvex airfoils at angles of attack near zero where the relax- 
ation time of the lower surface is close to that of the upper surface. 

Finally, the general case of an airfoil with different values of relaxation 
times on the upper and lower surfaces is discussed, and explicit expressions 
for the aerodynamic coefficients are given for the case where the condition 
on the upper surface is frozen (with the flow on the lower surface in nonequi- 
librium condition in general) . 


PRINCIPAL SYMBOLS 


a e 

a f 

c(£) 

c .p. 

CD 


■mi 


Ov 


equilibrium sonic speed 

frozen sonic speed 

a function defined by equation (4) 

Cm 

center of pressure, 77- 

C L 


drag coefficient, 
lift coefficient, 


drag 


(l/2)p co U 00 2 chord length 


lift 


(l/2) p^U^ 2 chord length 


moment coefficient, 


pressure coefficient, 


■Mlil - 3f2) 


pitching moment 

(l/2)p co U co 2 chord length* 

P ' Pco 


(1/2) PoA 


dh ) + (oh j (oq* 

^ P^oo \^P /oo J 


- 1 


2 



h 

lo 

K 


enthalpy 

modified Bessel function of the first kind, zero order 



c hord length 
K 


p pressure 

q nonequilibrium variable (e.g., degree of dissociation) 

s Laplace transform variable 

u,v velocity in the x and y directions, respectively 

free- stream velocity 

x,y Cartesian coordinates parallel and normal to the free- stream direc- 

tion, respectively 

a angle of attack (with respect of the free- stream direction) 



&(§) angle of the boundary at rj = 0 with respect to the free- stream 

direction 


M 


skew coordinates; £ 



and 



0 a function of £ = +0 U = 

p density 

q ratio of the relaxation lengths of the lower and the upper surfaces 
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1/2 


a 


a + s 
1 + s 


T o 

9 


0 


relaxation time behind the leading-edge shock 

normalized perturbation velocity potential, defined by the 

u - LL v u - ILj 

relations Kq> x = — — and K9y = — , or Cpg = — and 

U oo U oo U 00 


‘V' p l 


PfHo 


Laplace transform of cp> f°° exp(-s|)9( 


Superscripts and Subscripts 


( )* 
n 
( )’ 

( )u 

( h 
( >« 

( )nf 

( ) x , ( ) y n 
( ) 6 , ( ) t , J 


equilibrium value 

Laplace transformed quantity 

derivative of ( ) with respect to its argument 

upper surface of the airfoil 

lower surface of the airfoil 

free- stream value 

due to nonequilibrium free stream 

partial derivatives with respect to the subscripts 


PHYSICAL DESCRIPTION 


We shall consider an airfoil in a supersonic stream (fig- l) . The 
ambient stream is uniform but can be in a nonequilibrium condition. This can 
be expounded as follows : Let the state of the flowing gas be determined by 

the pressure p, density p, and a nonequilibrium variable q (e-g., degree 
of dissociation or vibrational energy). The enthalpy h, for example, is a 
function of these independent thermodynamic variables as h = h(p,p,q). The 
state of the gas is said to be in equilibrium if q satisfies the relation 
q = q*(p, p), otherwise the state of the gas is out of equilibrium. If q 
takes a fixed value, that is, q equals a constant, the state of the gas is 
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said to be in a frozen condition* In 
the problem we are considering, q is 
out of equilibrium everywhere including 
in the free stream, where 


± CL*(p ,P ) = 
00 00 




* 


The ambient flow is taken to be paral- 
lel to the x axis. It is assumed 
that the flow ahead of the airfoil is 
uniform with respect to both x and y. 
Thus , the relaxation time in the free 
stream, t^, is assumed to be large . 



Figure 1.- Airfoils in a nonequilib- 
rium free stream. 


Behind the leading- edge shock, however, 
the relaxation time 1 of the fluid is 
reduced due to compression, and chemical reaction or vibrational relaxation 
will take place. For flow past an airfoil at an angle of attack, the amount 
of reduction of the relaxation time across the leading-edge shock on the wind- 
ward surface is larger than that on the leeward surface because of higher 
compression on the lower side. We can classify the situation into three 
cases : 


1. When the angle of attack a is much smaller than the leading- edge 
slope 0 O of the airfoil, the relaxation time on the upper (leeward) surface, 
t u , is approximately the same as that on the lower (windward) surface, t^. 

The deviation from equilibrium is the same on the two surfaces. 

2. When a is smaller but of the same order as 0 O , compressive shocks 
still occur at the leading edge on both surfaces. Since < t u , however, 
the flow on the lower surface will be closer to equilibrium than that on the 
upper surface . 2 

3- When a is greater than 0 O , a compressive shock occurs at the lead- 
ing edge on the lower surface, but the flow on the upper surface is under 
expansion, which means t u is larger than hence, is very large. Conse- 

quently, the flow on the upper surface is in a frozen condition. 

We shall treat case 1 (t u = t^) in detail, giving both explicit expres- 
sions and numerical results. Only analytical results will be presented for 
the other cases, which will be discussed qualitatively. 

To attack this problem analytically, we shall restrict our study to a 
two-dimensional flow field that is perturbed slightly from the uniform 


1 The relaxation time is a characteristic reaction time. It is essen- 
tially the time required for a gas to reach an equilibrium state. 

2 An exceptional case is the flow past a wedge. The degree of nonequilib- 
rium in the free stream may be such that equilibrium on the upper surface is 
reached immediately behind the leading-edge shock, whereas it is not reached 
on the lower surface. This case is not treated here, however, since we are 
interested mostly in airfoils having closed profiles. 


5 



none quilibr ium free- stream condition. The relaxation time, however, is 
assumed to he reduced by an order of magnitude across the initial shock. Such 
an assumption for the relaxation time is hardly consistent with the assumption 
of small disturbances. When the disturbances are indeed small the relaxation 
time behind the leading-edge shock can be finite only if the free- stream 
relaxation time is finite so that nonuniform conditions exist in the free 
stream. When the disturbances are so large that the relaxation time is 
reduced by an order of magnitude, the assumption of small disturbances is vio- 
lated. Nevertheless, it is expected that for either of these two cases the 
linear theory should give indications of the flow phenomena on the airfoil 
surface, at least qualitatively correct. The nonequilibrium variable can be 
any of the common ones, for example, degree of dissociation, but we assume 
that only one nonequilibrium mechanism is of importance at a time. This is 
the same flow model with which Vincent! (ref. l) studied the flow over a wedge 
in a nonequilibrium oncoming stream. 

AN ANALYTICAL SOLUTION 


With the assumptions enumerated in the previous section we can describe 
the flow behind the leading-edge shock by a linear equation in terms of a 
perturbation velocity potential cp. This equation has been derived by 
Vincenti (ref. 1) for the upper plane, that is, for positive y. We shall 
obtain the pressure distribution on the upper plane, noting that the resulting 
expression is valid for the lower plane by simply replacing y by -y. In a 
linearized theory, the shock wave coincides with the initial Mach line. In 
the nonequilibrium linearized theory, the shock wave coincides with the ini- 
tial frozen Mach line. Since y = x/ is the initial frozen Mach line on 
the upper plane, all disturbances occur in x > P^y. It will be convenient 
to use, in place of the Cartesian coordinates x and y, the skew coordinates 
| = (x - P^y)/K and T] = P^y/K. The flow quantity K is proportional to the 
relaxation length TqU^ and is an important characteristic length in a non- 
equilibrium flow. To solve the resulting potential equation we shall use the 
method of Laplace transforms . 

The Laplace transformation of the flow equation is given by Vincenti, 
who obtained the transformed solution for^the upper plane, that is, positive 
T, (eq* (5*0 of ref. l) as 

4(s,tj) = L[cp(i,ri)] = A(s)e‘ (cr " l)sT1 + B(s)e (a+l)sT1 - ~ ^ (l) 

S P-£> 

The quantities A(s) and B(s) are to be determined from the boun dar y 
conditions. The quantity G^q^ - q^*)/(s 2 3 f ) comes from the nonequilibrium 
ambient condition. It vanishes when 00 = q^*. . Since and (q - q *) 
are normally positive (ref. l), the value of G^q^ - <1 *)/Pf is 5 nor mall y 
positive . 

The boundary conditions are: (l) the disturbances are bounded., and ( 2 ) 

the stream line is parallel to the boundary at the airfoil surface, which can 
be taken to be at y = 0 (hence r\ = 0) for the small disturbance theory. 



The first boundary condition means that the value of B(s) must be identically 
zero; otherwise the disturbances grow without bound as rj TO . The second 
boundary condition means 

9r|(i,0) - <Pg(i,0) - 8(|)/p f 


® n (s,0) - s«(b, 0) = &(s)/3 f - L[8(6)]/3f (£ 

where &(£) is the slope of the boundary at rj = 0 with respect to the free 
stream. This boundary condition gives the value of A(s). As a result we 
obtain the transformed value of cpt(i,Ti) = (u - U OT ) /U^ as 


L[q> 6 (!,Ti)] = s<d(s,ti) = 


G co^ - %>*) 1 _ _L 

Sff PfCT 


.(CT-X)TI _ goAo - Vj 

s3f 


The inversion of this transformed relation will give the value of 
(u - U OT ) /l^, which also gives the value of the pressure and enthalpy since, 
within the first order , 


Cp = 2(p - PJ/PJJ/ = -2(u - Uj/U* 


h - hoc = ~(u - UjU M 

Similar expressions for other flow quantities can also be obtained. 

The inverse transform of the last term in the right-hand side of equa- 
tion ( 3 ) is simply “G 00 (q co - q oo *)/Pf • No explicit inversion for the remaining 
terms in the right-hand side is known. Possibly they can be obtained by 
numerical means, (e.g., a method proposed by R. Bellman, et al. (ref. 2)). 

Our present interest is, however, only the pressure distribution on the 
boundary, that is, at rj = 0. For this case the exponential factor is -unity. 
The inverse transforms of l/(scr) and &/a can be obtained from a transform 
table (e.g., ref. 3)- They are 


-i 1 


c(|) = exp 


a + 1 


exp (" a 2 " ~ Io 


a - 1 


t dt 


where I 0 is the modified Bessel function of the first kind, zero order, and 


L' | = 8(0)c(|) + 


6* (t)c( | - t )dt 
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The pressure distribution on a boundary at t) = 0 with slope 5(|) is, 
therefore, given by the expression 

3f — ^ 29 =B(0)c(6)+ f 5»(t)c(| - t)dt + G (q - q *)[1 - c(g)] 

Poo U oo J C oo oo oo 

"which is quite general for any smooth &(£)• 

The above equation is derived for y (hence r]) = 0+. Therefore it gives 
the pressure on the upper surface. It is applicable on the lower surface 
(q = 0-) if one replaces y by -y in the derivation. This would result in 
replacing S by -8 in the above equation for the lower surface. Hence we can 
express the pressure on the upper surface, p u , and lower surface, p^, by the 
relations 

Pf — ■ - 2 °° = 5 u (0)cU) + f S u '(t)cU - t )dt + G (q^ - q *)[l - c(|)] 

oo oo O 

(6a) 

Pf — I s2 = -5 z (0)c(|) - r 5 z '(t)c(i - t )dt + G^q^ - q oo *)[l ~ c(|)] 

P U J o 


The first two terms in the right-hand side of equations (6a) and (6b) 
give the pressure distribution on the boundary when the ambient stream ahead 
of the leading-edge shock is in an equilibrium condition. The last term gives 
the effects of nonequilibrium free- stream conditions on the pressure. Within 
the first-order theory, therefore, the effect of nonequilibrium ambient condi- 
tion is in the form of a simple additive quantity. The consequences of such 
weak coupling to the aerodynamic characteristics of an airfoil will be dis- 
cussed in the next section. Because and q^ - q^ are normally positive, 
when the ambient stream is not in equilibrium, the pressure on the boundary 
is normally higher with the additive pressure proportional to q^*) • 

CHARACTERISTICS OF AIRFOILS FOR THE CASE K u = K z 


We shall now consider an airfoil at an angle of attack. First, we shall 
study thin airfoils in general. Later, we shall apply the results to a spe- 
cific airfoil, namely, the biconvex airfoil. The case for which the relaxa- 
tion times are approximately the same on the upper and lower surfaces will be 
treated presently. This is case 1 mentioned in the previous section. The 
case of K u ^ will be treated in the next section. 

The slope for the upper surface of an airfoil can be written, for 

convenience, 

e u = e(i) = e 0 f(i/i) 
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where Qq is the leading-edge slope, f ( |/z ) represents the variation along 
the chord line, and l is the chord length divided by K. (Recall that | 
is the streamwise distance from the leading edge divided by K. ) At an angle 
of attack the slope of the upper surface with respect to the free- stream 
direction is 6 U = 9 - a. The slope of the lower surface with respect to the 
free-stream direction is, for a symmetric airfoil, 8^ = -(0 + a). By substi- 
tution of 5 U an<3- into equations (6) we obtain the pressure coefficients 
for the upper surface, Cp u , and that for the lower surface, Cp_, , as 


gf c 

29 0 Cp n 


2 ^ 


0 O 


[1 - c(t)] + 


c(i) ,1 




c ( i - t ) dt 



(7) 


where the upper sign goes with Cp u and the lower with Cp^ . 

The various terms in equation (7) can he interpreted as follows: If the 

ambient stream is not in equilibrium (e.g., frozen in a certain degree of 
dissociation) , then after the flow passes the leading-edge shocks the sudden 
reduction of the relaxation time (due to the compression) induces recombina- 
tion , and the degree of dissociation tends toward an equilibrium value at a 
speed governed by the magnitude of the relaxation length behind the shock. 

This pressure rise induced by the recombination is given by the term inside 
the first bracket. It can also be, of course, due to other types of relax- 
ation, for example, molecular vibration. The term inside the second bracket 
gives the pressure distribution on the surface due to the contour of the air- 
foil. This portion of the pressure distribution corresponds to that of the 
"thickness case" (ref. 4) in classical aerodynamics. Here, of course, in con- 
trast to classical flows, the local conditions can be either equilibrium or 
nonequilibrium. Finally, the last term gives the pressure distribution due to 
the angle of attack, and corresponds to the pressure rise on a flat plate at 
an angle of attack in the classical case. With the expressions for Cp u and 
Cp^ we can compute the aerodynamic coefficients by the following relations: 
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while the center of pressure is 


c.p. = Cm/CL 


( 11 ) 


Using equation (7) for Cp u and Cp^ we obtain 

h Ct = i^o (s)di 


( 12 ) 


^rC m = 4 f ^U)di 


a 


(13) 


c.p. = — — 


f l ic(|)d| 


/•I 

i J c(i)di 


(14) 


4 C D'T 


t) c(s) + AiJ* (i) c( 5 ) U d5 


(15) 


where 


*(l) = c(|) j f'^JcU - t)dt 


Equations (12), (13)* and (l4) show that, within the first-order accu- 
racy, the lift, pitching moment, and the center of pressure of an airfoil are 
not affected hy the nonequilibrium condition in the ambient stream. They are, 
however, affected by the nonequilibrium condition behind the leading-edge 
shock. The explanation of such rather unexpected phenomena is that the pres- 
sure induced by the nonequilibrium ambient condition is the same on the upper 
and lower surface, hence, there is no effect on the lift and pitching moment 
(and, therefore, center of pressure). Such is not the case, however, for the 
drag . 


The drag coefficient consists of three parts: The first part is the same 

as a flat plate at an angle of attack and is given by the expression 


if c 

q q 2 Eflat plate 



c(£)dS 


'o 


( 15 a) 



The second part is the drag due the thickness of the airfoil and is 
given by the expression 



"^thickne 


ss 


4 

l 




c(| - t)dt 


d£ 


(15b) 


Finally, the effects of the nonequilibrium condition in the ambient stream to 
the drag is given by 


h 


°Dnf 



o(i)]di 


(15c) 


The lift-drag ratio can be obtained, of course, by use of equations (12) 
and (l5)- In particular, for the flat plate in an equilibrium free stream, 
the lift-drag ratio is simply equal to l/a. Thus, -within first-order accu- 
racy, the lift-drag ratio of a flat plate is not affected by the local non- 
equilibrium conditions - 

The values of Cp and Cp depend on the shape of the airfoil. To obtain 
representative values for them we shall consider a specific airfoil, namely, 
the biconvex airfoil, which is given by 

f = 1 - 2 


The corresponding Cp and Cp are 

Pf 

; ^ 1 


2 e 0 Cp u 


pf c 

200 Cp Z 


[1 - c(i)] + 


c(i) - -f f c(| - t)dt 


T f c(l) 

y o 


(16) 


h Cp = f f l |fe) + i 1 - 


( ; o i 


1 - 2 | ) [1 - c(i)n di 


(17) 


where 


^(i) = c(|) - | J c(£ - t)dt 
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Figure 2 shows the variation of lift, pitching moment, and (for a 
biconvex airfoil in an equilibrium ambient stream) components of the drag 
coefficients versus Z for a = 1.5* When the chord length is small compared 
with the relaxation length ( Z -* 0), the flow past the initial shock, not hav- 
ing sufficient time for relaxation toward the equilibrium condition, is 
essentially frozen. When the chord length is large compared with the relaxa- 
tion length (Z ■* °°) , however, the flow behind the leading edge is mostly in 
an equilibrium condition. Direct computation will also verify that in the 
two limits, the coefficients take the following asymptotic values: 



Figure 2.- Drag, lift, and pitching 
moment of a biconvex airfoil in 
a nonequilibrium flow field for 
(Mg 2 - l)/(M f 2 - 1) = 1.5. 

03 - 



Figure 3*- Frag on a “biconvex airfoil 
due to the nonequilibrium free stream 
for (Me 2 - l)/(M f 2 - l) = 1 . 5 . 


equilibrium: 


lim Cp 
l-+*> 


_4 _a 

JZ Pf' 


lim C m 

Z->-co 


_2 _a_ 

Pf 


lim Cp 

Z-* 00 



frozen: 


lim C L = ft, 
Z-*o 


lim C m 
1~*o 


2 a 

Pf 


li^Cp = A (f ♦ A 

Thus the coefficients in the frozen 
limit are larger than the correspond- 
ing ones in the equilibrium limit by 
the factor of 's/a. For the case of 
a = 1.5 ^ which is not unusual in mod- 
ern high- temperature facilities, this 
means a difference of roughly 20 per- 
cent . 


Figure 3 presents the change of 
Cp (ACD n f ) of a biconvex airfoil when 
the condition in the ambient stream 
is out of equilibrium. 

The drag induced by the nonequi- 
librium free stream is negative . The 
reason is that the pressure induced by 
the nonequilibrium ambient stream 
increases monotonically after the ini- 
tial jump at the leading-edge shock; 
consequently, a net thrust results. 
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Note that ACp^ is zero when the flow is frozen behind the leading edge 
since no recombination (for dissociating case) occurs. "When the condition 
behind the leading edge is equilibrium again zero because the 

pressure induced by the nonequilibrium stream is constant in the equilibrium 
limit . 

We see that, within the first-order perturbation approximation, the 
ambient stream nonequilibrium actually has no effect on the lift and the 
pitching moment. The local stream nonequilibrium, however, does affect the 
lift and moment, as well as the drag coefficients. The drag coefficient, on 
the other hand, is affected by the nonequilibrium condition in the free 
stream by an amount of 


“rtf = 


*“>0 °,x - 

P f l 




[1 - c(i)]d£ 


which is (l) proportional to G (q - q *) , (2) dependent on the shape of the 

y \ 00 00 CO 

airfoil, and(3J not dependent on the angle of attack. 


CASES WHERE K u ^ K z 


So far we have treated the case where the relaxation length on the upper 
surface of the airfoil is nearly the same as that in the lower surface (i.e., 
K u = = K, which is approximately correct when 0 o » a) . In general, due 

to the higher compression, the relaxation length is shorter on the lower sur- 
face, that is, < K u . If the angle of attack is larger than the leading- 
edge slope, K u > due to expansion. Since is assumed to be infinite, 

for this case, the condition is frozen on the upper surface. 

The case where Ku f is more complicated than the case we have 

treated. We shall briefly describe such a more general case, however, and 
conclude by giving the explicit relations for the aerodynamic coefficients 
when Qq < a, that is, K u 00 . 

For convenience, let 

K Z = K , and K u = ^ ^ K 

then 

= l , and Z u = pZ 

Note that 0 < p < 1. 

The relation in equation (6b) is still applicable directly for the lower 
surface since 6^ = 6 = x/K (for r] = 0) . For the upper surface the variable 
| must be replaced by = p£ and l by Z u = pZ . Thus we have 
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Pf 


Pf 


p 7 (6) - p 

L go __ 


P U 

H CO 00 


p u( g) - 


-5j(0)c(i) - / 5 z '(t)c(£ - t)dt + G (q^ - q. 00 *)[l - c(i)] 

J o 

(l8a) 


8u(0)c(|ig) + / 6 u ’(t)c(ni - t)dt + G (q - <^*)[l - c(n£)] 

( 18 b) 


where 


6u(n£) = 0(ni) - a = 0 o f ( )- 


a 


-5 Z ( £) = 0(|) + a = 0 o f( f )+ a 


(19a) 

(19b) 


The expressions for Cp are, therefore, 

|3~ G (q - q *) 

Cp U) = <S — [1 - c(nl)] + c(n£) 

dV Q U- Uq 

+ A f S f ’ ( A) 0(1X5 ' t>dt ■ fo C(tl4) (20a) 

^ %(|) = G “ ( \ V> f 1 - + <=< g > * T f f’(!)o( g -t)4t + £ 0(1) 

(20b) 


By direct substitution of Cp from equations (20) into equations (8) 
through (ll), one can obtain the aerodynamic characteristics of an airfoil 
when K u ^ . An interesting case is where \i = 0. For this case the flow 

on the upper surface is in a frozen condition- The value of Cp u is given 
by, therefore. 


^<W g > =-t + f 0) (21s) 

(Note that f < 1.) While Cp^ remains the same as before, that is, 

as % - v> t 1 - + t jT‘ Kt) o(i - *)« + f o o(6) 

(21b) 
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The corresponding aerodynamic coefficients are 


^ CL = f [1 . c(l)]+ ^ (s) . f fi 

' J Q i 


+ [c(£) + 1]| d| 
(22a) 


P f p 2 

T Cm T ,j 


Pf 


r a J%o - 


l a 



[i- c(|)] + -£ 


tU) - f ( 7 


20 o 1 rl 


3 C B=^C l+ ^ f^jc(|)-l + -° 

Uj x 


t(S) + f( | 


G (q. - q *) 

00 00 00 

a 


[1 - c(i)]Ui 


where 


+ [c(i) + 1]J> d§ 

(22b) 


(22c) 


tU) = c(6) + y I f \j) °(£ " t ^ dt 


( 23 ) 


It is apparent, from equations (22), that in the present case, because 
of the difference “between the relaxation lengths of the upper and lower sur- 
faces (K u ^ K i) , the aerodynamic coefficients are further affected. Consider 
Cp of equation (22a) as an example: Equation (22a) can he written as a sum 

of three terms , that is. 


Pp 

a 


C L 


P G (q - q *) 

Jr oo oo oo 

Z a 



[i - c(i)]di 


+ 


20, 


a 


o 1 
Z 




+ 


2 

Z 
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(24) 


Because the relaxation length is shorter on the lower surface (K^ < K^) , 
the pressure rise due to recombination (for dissociation equilibrium) is 
higher on the lower surface than that on the upper surface; the lift corre- 
sponding to this is given by the first term. 

The local flow condition is also different on the lower and the upper 
surfaces. The condition on the lower surface is closer to equilibrium. This 
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gives a higher pressure due to the thickness of the airfoil on the upper 
surface and a corresponding negative lift indicated in the second term. 

Finally, the suction on the upper surface due to the angle of attack is 
higher since the condition on the upper surface is closer to frozen. This 
gives an increase of lift due to the angle of attack as indicated in the 
third term in equation (24) . 


CONCLUDING REMARKS 


We may now make a few conjectures on the characteristics of lifting 
bodies in general, based on the observation of the preceding results. 

When the relaxation length is relatively constant on the "windward and 
the leeward sides, the effect of the nonequilibrium ambient condition on the 
lift and pitching moment is relatively small but will affect the drag and the 
pressure distribution. The local nonequilibrium conditions, however, will 
affect all of the aerodynamic characteristics. In fact, the aerodynamic coef- 
ficients are normally higher for bodies when the local stream is out of 
equilibrium. 

In the more general case the relaxation length is shorter on the wind- 
ward side than that on the leeward side. Thus, the condition on the windward 
side is closer to equilibrium, or, in other words, the leeward side is closer 
to frozen. In this case all of the aerodynamic coefficients are affected by 
both the ambient and the local nonequilibrium conditions. The ambient non- 
equilibrium condition will give an increase in lift. 
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